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Abstract. Horndeski models with a de Sitter critical point for any kind of material 
content may provide a mechanism to alleviate the cosmological constant problem. 
Moreover, they could allow us to understand the current accelerated expansion of 
the universe as the result of the dynamical approach to the critical point when it is 
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scenario that results when considering radiation and matter content, and conclude that 
their background dynamics is compatible with the latest observational data. 
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1 Introduction 

The cosmological constant problem [1, 2] is one of the main open problems in theoretical 
physics. This problem arises not only due to the ftne tuning necessary to (partially 
or completely) cancel the value of the vacuum energy, but also due to the need of 
repeating such a ftne tuning whenever a phase transition occurs changing the value of 
this vacuum energy by a ftnite amount or due to radiative corrections (see, for example, 
reference [3]). Moreover, the cosmological constant problem is present in any theory of 
gravity in which vacuum energy gravitates. One can think of a number of possible ways 
of cancelling the contribution of the vacuum energy in the dynamics of the universe, 
however, as Weinberg realised, this is a non-trivial task [1]. 

One of the possible stances towards this problem consists of admitting the ex¬ 
istence of some unknown symmetry that ftxes the value of the vacuum energy per¬ 
manently to zero, and consider alternative theories of gravity to describe the current 
acceleration of the Universe. The simplest extension of general relativity consist in con¬ 
sidering a four-dimensional metric theory of gravity with a scalar held non-minimally 
coupled to the metric. The form of the Lagrangian must be constrained in order to 
avoid the Ostrogradski instability, thus, one can consider a Lagrangian that only de¬ 
pends on first derivatives of the held. Galileon fields, however, are an example of a 
held Lagrangian containing second order derivatives which lead to a second order held 
equation [4]. The covariantization of this Lagrangian in curved spacetime with a held 
which is shift-symmetric was, therefore, a promising model [5]. However, this covariant 
galileons are strongly constrained by gravitational probes [6-8]. 

Generalized galileons [9] have evolved beyond the covariantization of reference 
[5] and comprise models with the most general Lagrangian containing second order 
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derivatives of the field and leading to second order field eqnations, withont assnming 
any symmetry for the field. The resnlting Lagrangian has been proven [10] to be 
eqnivalent to the already known (although also slightly forgotten) Horndeski theory 
[11]. In the context of this theory one could look for an adjustment mechanism able 
to alleviate the cosmological constant problem. The fab four models [12, 13] avoided 
Weinberg’s no go theorem by relaxing one of its assumptions. The resulting models 
have a scalar held which is able to self-tune, that is, the held screens the contribution 
of any value of the vacuum energy leading to Minkowski. Considering that such a 
self-tuning has to be dynamical to allow non-trivial cosmological dynamics, one is led 
to the conclusion that the fab four cosmologies have a Minkowski critical point for any 
value of the vacuum energy and any kind of matter [14]. 

In reference [15], the concept of self-tuning was extended to consider that the 
held screens the vacuum energy leading to a de Sitter instead of a Minkowski vacuum. 
Thus, in this scenario, the current accelerated expansion of the universe results from 
the dynamical approach to a stable critical point. It has been emphasized that the 
value of A characterizing the critical point is completely independent of the vacuum 
energy and, therefore, its value is unahected by phase transitions. However, in order 
to alleviate the cosmological constant problem it would be necessary that not only the 
critical point is an attractor but also that a long enough radiation and matter phases 
can be described before the attractor is approached, even if the vacuum energy takes 
a huge value. 

As found in reference [15] there are two families of models able to self-tune to a 
spatially flat de Sitter vacuum. The first family has a minisuperspace Lagrangian den¬ 
sity that is proportional to the derivative of the held. The cosmological consequences 
of these linear models have been studied in reference [16]. Although promising re¬ 
sults have been presented, the considered models were not completely able to describe 
the background cosmology of our Universe, nor results concerning the stability of the 
critical point can be obtained in the general linear case. 

The second family of models has a minisuperspace Lagrangian with a non-linear 
dependence on the temporal derivative of the held which vanishes when evaluated at 
the Hubble parameter characterizing the critical point. These models contain three 
arbitrary functions of the held and its derivatives. Focusing our attention on shift sym¬ 
metric functions can simplify the general treatment of this family of models. Moreover, 
it has been proven that the galileon interactions are nonrenormalizable and that the 
additional terms generated by quantum corrections contributing to the Lagrangian are 
strongly suppressed at energies below the cut-oh [4, 17, 18]. Therefore, we expect that 
the screening mechanism, which is a result of the particular form of the Lagrangian 
[15], would not be spoiled by held and/or matter quantum corrections for the shift- 
symmetric models (which are curved spacetime generalizations of the galileon). It must 
be noted that by assuming a shift-symmetric held, one is not automatically in the co¬ 
variant galileon case ruled out by observations [6-8] and, therefore, these models are 
interesting and worthy of being explored. In this article we study the cosmology of the 
shift-symmetric non-linear Horndeski models with a de Sitter critical point regardless 
of the material content. 
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This article can be outlined as follows: In section 2 we summarize some character¬ 
istics of the Horndeski cosmological models able to screen any value of the cosmological 
constant to a given de Sitter space, focusing attention on the case of a shift-symmetric 
held. In section 3 we study the dynamical system for these models showing that the de 
Sitter critical point is always an attractor if the material content of the universe satis- 
hes the null energy condition. In section 4 we consider some simple cases which allow 
us to gain some intuition regarding the background cosmology of these models. This 
study allows us to consider a particular satisfactory model which we present in section 
5. In section 6 we summarize the main conclusions of the article. As the explicit form 
of the general Lagrangian of the satisfactory models can be of great interest in order 
to consider further physical scenarios, we include in appendix A and B the functions 
which appear in this Lagrangian according to the notation of reference [11] and [9], 
respectively. 


2 The non-linear models 


The minisuperspace Lagrangian of the Horndeski models with a spatially flat de Sitter 
critical point for any value of the vacuum energy can be expressed as [15] 

L = Teh + L linear “1“ -^nl H” -^m* (2.1) 


In this Lagrangian we have explicitly written the Einstein-Hilbert term 

Teh ~ 

which belongs to the hrst family of models^ 


.^linear ® E [37X17,(0)+ 0 IE,(0) 




i=0 


The second family of models has a Lagrangian given by 


3 

= ( 0 , 0 ) H\ 

i=0 


( 2 . 2 ) 


(2.3) 


(2.4) 


The functions Ui^cf)), 117(0), and Xi (^0, 0j are arbitrary up to the following constraints 
which ensure the existence of a de Sitter critical point characterized by A, the late-time 
value of the Hubble rate, for any kind of material content [15]. These are: 

3 3 

(2.6) 

i=0 j=0 

^It must be noted the E-H term is not able to self-tune by itself [15]. 
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and 


3 

0 ) a */2 = O . ( 2 . 6 ) 

i=0 

The first constraint ensures that the Lagrangian density evaluated at the critical point 
has the form required to allow the held to self-tune, which is 

Cep = 3a/A h{(f)) + 0 (2.7) 

with h(0) an arbitrary function. The second constraint forces any non-linear depen¬ 
dence of the Lagrangian on 0 to vanish at the critical point to agree with equation 
(2.7). Furthermore, we consider that the material content is minimally coupled to the 
metric and non-interacting with the held. Thus, we have 

Am = -a^p{a) = -a^'^ps{a), with p^(a) = ( 2 . 8 ) 


The cosmology of the linear models was considered in reference [16]. In the present 
article we instead study the cosmology of the non-linear models, that is, we assume 
that Ui{(j)) = lTj(0) = 0. The held equation for the resulting models is given by [15] 


E 


A'i,^ - 3A'.. iff - lA, - A, ij - A, 




JJ 




w = 0 . 


(2.9) 


where one needs at least one i such that ^ 7 ^ 0 with i 7 ^ 0 to obtain a dynamical 
self-tuning, that is, a term depending on H in the held equation. Moreover, taking 
into account condition ( 2 . 6 ) and its derivatives, it can be checked that this equation 
trivially vanishes at the critical point, that is, when H —)■ a/A. On the other hand, the 
Hamiltonian density is 

H = HEH + 'Hnl + Hm = 0, (2.10) 

where 


77eh 

•Hnl 

Hm 


3 

5^[(i-l)A.(.^, (^,0 

j=0 

P(a), 



( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


respectively, and we have dehned for obvious reasons. Since using (2.6) and its 
derivatives gives 

3 

= ( 0 , 0)aA2, (2.14) 

i=0 

and we are assuming at least one ^ 7 ^ 0 with i 7 ^ 0 , this implies that 

3 

(2-15) 

j=0 
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and consequently, 7 ^ 0 at the critical point. Thus, the field is able to self-tune 

[15]. The modified Friedmann equation is given by (2.10). It can be explicitly written 
as 

If we define = p/(3MpjTr^) and = p^/{3MpiH^), with given by equation 
(2.12), we can rewrite the modified Friedmann equation as 

n + n^ = i. (2.17) 

Finally, one can define an equation of state parameter associated with the field as 
= Pep!P 4 >-, which can be written as 

and an effective equation of state parameter as 

Weff = -1--^. (2.19) 

2.1 Shift-symmetric case 

Let us now focus our attention on models that are invariant under shift redefinitions 
of the field, 0 —)■ 0 -|- c. In this case the functions Xi are just a function of the field 
derivative, Xj(0). Let us redefine the functions as 

u.i’) = 14 (2'20) 

for later convenience. In terms of these new functions, conditions (2.6) and (2.15) are 
respectively given by 


3 3 

^/i( 0 ) = O, and 4 ,^( 0 ) ^ 0 . ( 2 . 21 ) 

2 = 0 2=0 

As the Lagrangian is invariant under 0 —)■ 0 -|- c, there is a conserved quantity Sq with 

^ 1 (97 

= = ( 2 . 22 ) 

i=0 

where h = H/\fK. If Sq = 0, taking into account condition (2.21), one can conclude 
that the universe stays either in a Minkowski state, h = 0, or in a de Sitter state 
with, h = 1 (as — 0); throughout the whole evolution. For the more general 

case with Sq 7 ^ 0 , one can obtain from equation ( 2 . 22 ) that either h —)■ 0 or h —)■ 1 
when a —)■ cxD. Thus, for these models a late time de Sitter evolution is necessarily 
characterized by h = 1 (77 = V^). 
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As the functions fi{(f>) only depend on 0, the held equation is greatly simplihed 
and can be expressed in terms of V’ = 0- Dehning now, N = lira, as the new time 
variable, and denoting with a prime the derivatives with respect to N, this equation 
can be written as 

3 3 3 

3h ^ + hip' ^ + h' ^ ifi^^{ip)h' = 0. (2.23) 

2 = 0 2=0 2 = 0 


The modihed Friedmann equation is just the constraint (2.17), which can be expressed 
as 


12 T — 1, 


(2.24) 


where 


3 


i=0 

Finally, as each specie is conserved, we have 


(2.25) 


12 = 12s 

s 


with 


12 ' = - 12 , 


. h' 
3(1 + Wg) + 2 — 


(2.26) 


Thus, the total density parameter satishes the following differential equation 


12 ' 


-12 


3(1 Tw) + 2 


h'' 

~h 


(2.27) 


where we have dehned the total equation of state parameter of matter huids as 


1 + tc 




(2.28) 


which should not be confused with Wes- This quantity, w, generically is a function of 
the scale factor whenever there is more than one specie. For the cases when Wg is a 
constant, equation (2.25) can be integrated to yield 


12 


S 


12s,* K 


exp [—3(1 + Wg){N 


N.)], 


where the subscript * means evaluation at the initial condition. 


(2.29) 


3 Dynamical system analysis for shift-symmetric cases 

Considering the equations presented in the previous section, it can be noted that we 
have three functions of N, those are ip, 12, and h, which appear with their hrst deriva¬ 
tives in the held equation (2.23), the Friedmann equation (2.24), and the conservation 
equation (2.27). Moreover, the Friedmann equation (2.24) does not contain derivatives 
of any of these quantities, thus it can be seen as a constraint which can be used to 
decouple the dependence on the derivatives of equation (2.23) and (2.27). 
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With this aim let us hrst rewrite the held equation (2.23) as 



^7'Pi (h, V’) + h'P 2 (h, ij) + Po (h, ij) = 0, 

(3.1) 

with 

3 

Po(h, 

i=0 

(3.2) 


3 

Pi (h, Pj) = 

i=0 

(3.3) 


3 

P 2 (h, ?/>) = 

i=0 

(3.4) 

Given that G = 1 —with given by equation (2.25), and substituting the resulting 

expression and its derivative in equation (2.27) we obtain 


ip'Qi (h. Ip) + h'Q2 (h, pj) + 3(1 -h w)Qo (h, ip) = 0 , 

(3.5) 

where 

3 

Qo {h, ^Z^) = -1 + ^ [(i - l)fi + ipfik = Gp - 1 

i=0 

(3.6) 


3 

Qi (h, i’) = Y1 + V'/i.yy] 

i=0 

(3.7) 


Q2{h, Ip) = h~^ |-2 + [(i - l)/i + . 

(3.8) 

From equation (3.1) and (3.5) one can obtain 



, 3{1 + w)QoPi — QiPq 

Q 1 P 2 — Q 2 P 1 

(3.9) 

Assuming a single-component universe with w constant, these two equations form an 
autonomous closed system. Taking into account ( 2 . 21 ), it can be noted that Pq (1, V’) = 
Pi (1, V’) = 0 and P 2 (1, V’) = Qi (1) V’) 7^ 0) cind, therefore, P = 0 and ip' = —3(1 -|- 
w)Qo/Qi at h = 1 . Thus, we have a de Sitter critical point for any material content 
(that is, any w) provided Qq (1, V’c) = 0 and Qi (1, -^c) 7^ 0. The critical point is, 
therefore, characterized by 


3 

{h = 1, Ip = ipc}, with fi{ipc) = 1, 

i=l 

(3.11) 
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obtained by demanding (5o(l, V’c) = 0 and using the Ihs relation of (2.21) and its first 
derivative, and 

3 

(3.12) 

i=l 

which results from requiring Qi (1, ^ 0 and using now the the second derivative of 

the Ihs of (2.21). Taking into account equations (2.24) and (3.6) it can be seen that 
this corresponds to flc = 0, as it should be expected. 

Let us now study the stability of the critical point using linear stability theory. 
The eigenvalues of the Jacobian matrix of the system given by equation (3.1) and (3.5) 
evaluated at the critical point (3.11) are 

Ai =-3, Aa =-3(1+ ■«;). (3.13) 

Thus, the critical point is an attractor if 1 Try > 0, a saddle point if 1 + ry < 0, and one 
should go beyond the linear stability analysis for pure vacuum. It must be emphasized 
that in this section we are considering a single component universe with constant w. 
Nevertheless, this result suggests the stability of more general models with components 
satisfying the null energy condition. 


4 Recovering early time GR with functions of the same order 


In order to understand which class of models can deliver an interesting cosmological 
behaviour, we assume in this section, that the functions fipip) either vanish or are of 
the same order at early times. In particular, we seek to identify which models possess 
a general relativistic cosmological history at times corresponding to very large values 
of the Hubble parameter, h ^ 1. One can note that, due to constraint (2.21), models 
with only two non-vanishing functions, these are automatically of the same order as 
we must have fi{'ip) = ~fjidP)- Taking into account equations (2.19) and (3.9), the 
effective equation of state parameter is 


l + WeS 


2 3(1 -|- w)QoPi — QiPq 
3h Q 1 P 2 — Q 2 P 1 


(4.1) 


Thus, the referred assumption allow us to study the regime h S> 1 by keeping only the 
terms which are higher order in h in the functions Qfs and Pfs. We can classify the 
different cases by considering the non-vanishing /j with the largest value of i. 


Case I. One can see that if the non-vanishing function with the largest value of i is 
fsi'f’), one can approximate equation (4.1) by 


. I ^ 2 (1 -|- w) ( 2/3 -|- ' 0 / 3 ,p) fs^tp^ — ( 3 / 3 ,p + ' 0 / 3 ,pp) /s.p 
3 ( 2/3 + '0/3,-i/>) f3,ipi> — ( 3 / 3 ,p + '0/3,pp) /3,p 


(4.2) 


for h S> 1. Thus, assuming 0 < tc < 1 one has two different extremal cases, either 


1 + Wes - ^{3 + w), for 


(2/3 + "0/3,p) /3,pp| S> I (3/3_.0 -I- /3,p|, (4.3) 




or 


2 

1 + Weff ^ g, for I ( 2/3 + < I ( 3 / 3 ,^ + V'/s.i/;^) / 3 ,^|- (4.4) 

Neither of these two cases is suitable to describe our Universe as radiation and later 
non-relativistic matter are expected to drive its dynamics at early times. One cannot 
obtain intermediate cases compatible with our cosmological history without any hne- 
tuning either. 

Case II. Let us now consider that the non-vanishing function with the largest value 
of i is Considering the regime h ^ 1 in equation (4.1), we have 


1 I ^ (1 -f w) (1 — /2 — '0/2,p) /2,yy + ( 2 / 2 ,.,/, + '0/2,yy) /2,y /. -n 

(1 — /2 — '0/2,y) /2,yV’ + ( 2 / 2 ,y + '0/2,yy) /2,y 
For 0 < tc < 1 the extreme regimes are 

1Weff ^ 1 +for |(l-/2-0/2,y)/2,yy| > |(2/2,p + '0/2,yy)/2,y|, (4.6) 


and 


1 + Weff~l, for |(l-/2-0/2,y)/2,yy| < |(2/2,p + '0/2,yy)/2,y|. (4.7) 


The case given by equation (4.6) corresponds to a cosmological evolution which could 
be compatible with a general relativistic behaviour at early times. Analysing this case 
in further detail, we require a positive and small value of for h ^ 1. Taking into 
account equation (2.25) one has that 


~ /2 + 0/2,V, > 0, 


(4.8) 


for h ~ h* 3> 1 and 0 ~ 0*. The smallness of this quantity can be controlled 
by imposing a given held initial condition, 0*, when integrating equations (3.9) and 
(3.10). From inequality (4.8) we must have 0 < 1 — (/2 + 0/2,1/;) < 1 at early times, 
therefore, condition (4.6) is satished if 

|/2,V'i/;| s> I (2/2,,/, + 0/2,i/;i/;) /2,i/;|, (4.9) 

for 0 ~ 0*. It is easier to ensure that conditions (4.8) and (4.9) are satished close 
to 0* if the held only takes either positive or negative values. A sufficient condition 
to guaranty this is to set a positive (negative) initial condition with 0 ' > 0 (< 0 ). 
Evaluating (3.10) for h S> 1, and taking into account condition (4.6), one has 

0' ~ 3w^^. (4.10) 

72,yp 

Thus, if0* > 0, one should require sign(/ 2 ,^) = sign(/ 2 ,^^), and sign(/ 2 ,p) 0 sign(/ 2 ,^p) 
otherwise. One has, of course, to check that if 0* > 0 then 0* < 0c, and 0* > 0c if 
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Figure 1. Illustration of a case II dynamics. Model given by the functions (4.11) with 
a = 10“^, 6 = 2, c = 2a and d = 26. Initial conditions have been taken at equivalence, those 
are z* = Zeq = 3.3 x 10^, 6,* = /igq = 1.6 x 10®, and = 0.00001. On the upper panel we 
show the evolution of the effective equation of state parameter for our model, depicting the 
current value favoured by observations. On the bottom panel we show the evolution of the 
0-parameters. 


"0* < 0, with 'ipc given by equation (3.11), to ensure that we are in the solution branch 
with the attractor. 

These models are, therefore, very promising. Nevertheless, by studying different 
cases with power functions or exponential functions, we have concluded that their late 
time behaviour is not completely satisfactory. In hgure 1 we show the result of the 
numerical integration of equations (3.9) and (3.10) for a model described by 

= = fq{^) = -ae^^-ce^'^, (4.11) 

which has a behaviour similar to a model with power functions. It must be noted 
that the parameters of the model have to be such that conditions (4.8) and (4.9) are 
satished, with b and d of the same order. Both kind of models present a current value 
of Wes too small to be compatible with observations, as it can be seen in hgure 1. One 
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Figure 2. Case II dynamics. Model given by the functions (4.11) with a = 9 x 10“^, 
b = 2, c = 2a and d = 2b. Initial conditions have been taken at equivalence, those are 
z* = Zeq = 3.3 X 10^, /i* = heq = 1.6 X 10®, and V’* = 0.00001. 

can only find a viable valne for taefr,o at the price of introdncing large amonnts of early 
dark energy. This is shown in fignre 2. 

Case III. The remaining case necessarily corresponds to the sitnation with only two 
potentials, /i('0) and /o('0)- In the first place, taking into acconnt eqnation (3.9) into 
eqnation (4.1) for h 3> 1, we get 


1 + Weff — 1 + tc, (4-12) 

for any form of the fnnctions, making this models of particnlar interest. In the second 
place, from eqnation (2.25), we mnst impose 

~ ^ > 0, (4.13) 

for h ~ h* 1 and -0 ~ '0*, where we have implicitly assnmed that /i/('0/i,p) 
evalnated at "0* is not arbitrarily large. The smallness of this qnantity can be tnned 
by fixing the valne of 0* once a particnlar form of the fnnction j] is given. Condition 
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(4.13) is immediately satisfied if sign('0) = sign(/i_^), for %Ij 'll )To ensure that ip 
does not change sign in the neighbourhood of V’* one can impose V’ > 0 with pj' > 0 
for h ^ 1 or alternatively ip < 0 with ip' < 0. This is, of course, a sufficient condition 
although not necessary. Considering the regime h 3> 1 in equation (3.10) one gets 


3(l-n;)/i,p ^ 

2/i,v>p 


(4.14) 


Thus, for 0 < w < 1, -0' > 0 (< 0) if sign(/i,^) ^ sign(/i,^^) (sign(/i,^) = sign(/i,^^)). 
Taking into account also condition (4.13) and the value of the field at the critical point 
(3.11), we can restrict our attention to consider models with either 0 < ip^ < ip^, 
/i,p > 0 and < 0, or ip^ < ip^ < 0, < 0 and < 0, with ip^ given by 

equation (3.11). 


w = 0 



h h 


Figure 3. Phase space diagram for case III. In the left panel we show the results for the 
model with a function given by equation (4.15) with c = 2, b = 3 and a = 1 and d = 2b. In 
the right panel we depict the model corresponding to equation (4.16) for b = a = 1/2. The 
dynamics evolves away from the de Sitter critical point (indicated by the filled circle) in the 
branch of solutions with large initial values of h. 

Despite the simplicity of these models, and encouraging characteristic shown in 
equation (4.12), the arguments suggesting an early time general relativistic behaviour 
cannot be taken for granted if unexpected features appear in the phase space. This is 
precisely what we hnd when studying particular models of this kind. Assuming ip > 0 
one can consider a model given by 

h{iP)=c-be-'^^, (4.15) 

and a second model with 

f,{ip) = br, (4.16) 
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being /q = — /i, c, b, a > 0. It must be noted that for the hrst case one needs 
l<c<l + 6to have 'ipc > 0. The bounds come from equation (3.11), since c = 
1 + < 1 + 6 for aipc > 0 and 6“^^= = 6/(c — 1) > 0 with 6 > 0. Looking 

at the phase diagram of these models in hgure 3, one can conclude that even if the 
critical point is an attractor, solutions with a large value for the initial condition for the 
Hubble parameter, h*, correspond to a branch for which -0 continues to grow. Thus, 
cosmologies with viable initial conditions will not reach the de Sitter attractor. 

5 Beyond the simplest assumption 

Under the assumption considered in the previous section, that is, that at early times 
the functions either vanish or are of the same order, we have been able to obtain 
promising results for models in which the function with larger value of i is / 2 (' 0 )- 
These models have interesting early time cosmology, even though the current value of 
the cosmological parameters differ from values suggested by the current observational 
data. In an attempt to obtain a more satisfactory cosmology at present time we will go 
beyond the assumption taken in the previous section still considering that / 2 ('?/’) is the 
non-vanishing functions with the largest value of i and that it satishes the conditions 
presented in the previous section. It must be noted that, because of condition (2.21), 
we need to have two additional functions to avoid them to be of the same order as 
/ 2 ('*/’)• We consider a model with only /2('0) and /i('0), with an extra term that 
modihes /i('0) such that it differs substantially from / 2 (' 0 ) at least in some range of 
V’. Thus, we study the following model: 

= ai)^, /i(?/>) = /o(V’) = (5.1) 

Looking to the past, if h increases faster than 'ip decreases, the terms '0"' dominate and 
this models has a consistent early cosmology for small enough values of a, as in case II 
of the previous section. If p) decreases faster than h increases back in time, however, the 
models could develop the same characteristic in the phase space as shown in hgure 3. 
(We have studied that this can happen for large values of a.) The functions appearing 
in the general Lagrangian (without restriction to the minisuperspace) of this model 
are included in appendix A and B, depending whether one considers the Horndeski 
formulation [11] or the Deffayet et al. expression [9]. 

Considering a universe hlled with radiation and non-relativistic matter, one can 
numerically integrate equations (3.9) and (3.10). Taking into account the evolution of 
the conserved matter content, equation (2.29), one can depict the relevant quantities 
of the model. The evolution of these models is quite satisfactory. For example, for a 
model with 

a = 10"^ n = 2, m = 1, /S = 7.5, (5.2) 

taking hgq = 1.6 x 10® and 'i/'eq = 1.0010074751 at equivalence Ze^ = 3.3 x 10^, we 
obtain consistent evolutions of the fl-parameters and current values compatible with 
observations, = 0.300929 and = 9.11907 X 10 ®, and small quantities of early 
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z 


Figure 4. Model given by the functions (5.1) with parameters (5.2). Although the evolution 
of other functions of the model is satisfactory, the decrease of the field equation of state 
parameter at the present time, characterized by Wa, seems too large to easily agree with 
the observational data. 

dark energy, fltp^rec = 6.77705 x 10“^. Nevertheless, considering a Taylor expansion of 
the held equation of state parameter, depicted in hgure 4, around the current scale 
factor as 

~ Wo + Wa(l - a), (5.3) 

we obtain wq = —0.945417 and Wa = 1.49429, which are at the very best only in 
marginal agreement with the observational constraints (see e.g. [19]). On the other 
hand, for a model with 

a = 10“®, n = 16, m = 1, 13 = 1.5, (5.4) 

the cosmological history becomes entirely compatible with current data, considering 
conditions heq = 1.6 x 10^ and V’eq = 1.03999 at equivalence Zeq = 3.3 x 10^. The 
evolution of the energy densities for this particular case is shown in hgure 5. The 
results are compatible with current observational data, as we obtain Om,o = 0.302601 
and Qr,o = 9.16972 x 10“®, avoiding early dark energy with fltp^rec = 3.90298 x 10“^. In 
hgure 6 we show the evolution of the ehective equation of state parameter and the held 
equation of state parameter. As we show, the current value of Wes is compatible with 
the data. Moreover, we obtain for our model Wq = —0.980623 and Wa = 0.0436749, 
equation (5.3). It can be verihed that this model presents a future brief phantom epoch 
before approaching the cosmological constant behaviour. 

6 Summary and further comments 

In this article we have consider Horndeski cosmological models that may alleviate the 
cosmological constant problem by screening any value of the vacuum energy given by 
the theory of particle physics. In particular, we have studied in detail the non-linear 


14 




-Iogio(l+z) 



Figure 5. Model given by the functions (5.1) with parameters (5.4). On the top panel we 
compare the evolution of the fl-parameters, and on the bottom panel, the matter and field 
content in terms of the logarithmic of the energy densities, where n = 87rG/(3A). 


family of models obtained in reference [15] which have a de Sitter critical point for any 
material content. Furthermore, we have considered models protected by a particular 
symmetry, a shift symmetry of the held. 

As we have shown, for these shift-symmetric models the de Sitter critical point 
is indeed an attractor. Thus, we can understand the current accelerated expansion of 
our Universe as the result of the dynamical approach of the held to the critical point, 
being the value of —)■ A at this point completely independent of the vacuum energy. 

The background cosmological evolution of the models studied in this article sug¬ 
gests that these models are in even better footing than the linear family considered in 
reference [16]. It is clear from our analysis that there is a region of parameter space 
that it is incompatible with current observational data, hence, these models are sus¬ 
ceptible of being ruled out. However, we have also identihed a particular case able 
to describe currently available observational data which depends of four parameters. 
Imposing the integration conditions at matter-radiation equality, this model provides 
us with a value of the density parameters and equation of state held parameter at 
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-log 10 ( 1 +z) 
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Figure 6. Model given by the functions (5.1) with parameters (5.4). On the top panel we 
show the evolution of the effective equation of state parameter and the current value favoured 
by observations. On the bottom panel we show the evolution of the field equation of state 
parameter. 


the present time within observational bounds. Moreover, the contribution of the held 
at early times is negligible, recovering a cosmological dynamics compatible with that 
produced by general relativity. In particular, there is no early dark energy. 

In order to scrutinise these models we are now required to face them against 
observables that depend on the evolution of the held and matter huid huctuations. 
Another possible extension of the current work consists in investigating how the linear 
and non-linear contributions to the minisuperspace Lagrangian ahect the dynamics 
when they are both present. This will be carried out in future work. 

As we stated in the introduction, these models can alleviate the cosmological 
constant problem only if a long enough radiation and matter phases can be described 
before the attractor is approached for any value of the vacuum energy. This point has 
not been addressed in the present work. One could expect the dynamical screening 
to start before complete screening has taken place at the critical point. Moreover, for 
some models it may even be possible that the screening could be more effective for 
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vacuum energy than for the material content, allowing a consistent cosmology. Such 
study should, however, be carefully carried out in a follow up project. 
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A Horndeski functions 

In this paper, we have started by including the Lagrangian already restricted to the 
minisuperspace, which is the case of interest for studying the background cosmology of 
a particular model. Nevertheless, once a particular satisfactory model has been found, 
one needs to write the general Lagrangian to study other consequences of the model. In 
order to facilitate that study for future works, we include here the general Lagrangian 
of the model presented in section 5. 

The Horndeski Lagrangian, as expressed in reference [11], can be written as 

o 

— 3 [2E^ -|- X Kg] -f- Kg, (A.l) 

where 

F,x = ^1,0 — — 2Xk3,x, A-2) 

X = and k^ (0, X ) are arbitrary functions. We are considering shift- 

symmetric models, therefore, the functions are just dependent on the kinetic term 
Ki { X ). As it was hrst shown in reference [13] for the minisuperspace Lagrangian 
in the general case, the functions of Lagrangian (2.4) are related with the functions 
appearing in Lagrangian (A.l) through 


Ag — 

(A.3) 

Xl = —3 Qy -I- ^ 7 ^ 00 , 

(A.4) 

X 2 = 12 FxX — 12 E, 

(A.5) 

A 3 = 8 Ki,x 0^, 

(A. 6 ) 


(A.7) 


£h = ys 
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where we have simplihed the expressions due to the shift-symmetry. Thus, the La- 
grangian for the models given by (5.1) have the following Horndeski functions 


— Cl, 

«3 = 




4(n - 2){n - 1) 


a ^ + C2 + 


C3 




1/2 ■ 


F = 


2{n - 2 ) 
Kg = Mp^y/X 




n 


n — 3 




m 


m + 3 




C 4 , 




(A.8) 

(A.9) 

(A.IO) 

(A.ll) 

(A.12) 


being the Cj’s integration constants. It must be noted that the terms appearing multi¬ 
plied by Cl, C 2 and C 4 can be combined in a total derivative; therefore, these constants 
can be hxed to zero. The constant Cg does not appear in the minisuperspace La- 
grangian, therefore, the corresponding term is not able to self-tune to de Sitter by 
itself although it does not spoilt screening (as it happens in the case of the linear 
models with two potentials [15]). 


B Deffayet et al. functions 

Deffayet et al. independently found the Horndeski Lagrangian in reference [9], ex¬ 
pressed in a form which is currently more used in the literature. Assuming a shift- 
symmetric field, this is 

C = K{X) - Gs{X) □</. + GXX)R-2 G^,x [(□</-)' - (V^V.c/.)'] 

+ G5(X)G^,V^V> + [(00)' - 3(n0)(V^V.0)2 + 2(V^V.0)='] . (B.l) 

Taking this symmetry into account, the dictionary hrst presented in reference [10] 
relating the former Lagrangian with Lagrangian (A.l) can be expressed as 


1 ^ 9 , 

(B.2) 



Xks — / dX'ftg, 

(B.3) 

2F 2Xk3, 

(B.4) 

— 4ki. 

(B.5) 
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For the model given by equation (5.1), taking into account equations (A.8-A.12), we 
have 


K 


-3Mji AI3 


(B.6) 

(B.7) 





m + 1 


(B.8) 

(B.9) 


where we do not take into account the terms leading to total derivatives, and we have 
dehned c = 2 C 3 . As c does not affect the background cosmology, the value of this 
constant is not restricted by our analysis, and it can be fixed to the more convenient 
value. If one considers the limit case a = 0, one is in case III of section (4) for m < 0. 
This case is particularly simple as the Lagrangian (B.l) only contains two terms. On 
the other hand, the limit case f3 = 0 corresponds to case II of section (4). 
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